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Abstract 

Our goal is to discover possible new 4-dimensional euclidean solutions (instantons) in 
fundamental SU (2) Yang-Mills-Higgs theory, with a constraint added to prevent collapse 
of the scale. We show that, most likely, there exists one particular new constrained 
instanton (I*) with vanishing Pontryagin index. This is based on a topological argument 
that involves the construction of a non-contractible loop of 4-dimensional configurations 
with a certain upperbound on the action, which we establish numerically. We expect I* to 
be the lowest action non-trivial solution in the vacuum sector of the theory. There also 
exists a related static, but unstable, solution, the new sphaleron S*. Possible applications 
of I* to the electroweak interactions include the asymptotics of perturbation theory and 
the high-energy behaviour of the total cross-section. 



1 Introduction 



We have conjectured JT| the existence of a new constrained instanton in the vacuum 
sector of euclidean SU(2) Yang-Mills-Higgs theory. Our argument consists of two 
steps : 

(1) the construction of a non-contractible loop of 4-dimensional configurations of 
the Yang-Mills and Higgs fields, starting and ending at the classical vacuum; 

(2) the proof of a certain upperbound on the action over this non-contractible loop, 
which is crucial for having a genuine new solution, as will be explained later on. 
For the last step, which is by far the most difficult of the two, we have to resort 
to numerical methods. The numerical results of our previous paper |]J were not 
entirely conclusive in establishing this upperbound on the actionprofile. Here, we 
present further numerical results that are conclusive, at least for certain values of 
the parameters in the theory. 

Compared to our earlier work we have made the following improvements. First, 
the constraint on the scale is treated dynamically, whereas before the scale was 
fixed by hand. Second, the ansatz of the non-contractible loop is generalized, in 
order to give the fields more freedom to relax to a lower value of the action. Third, 
the efficiency and accuracy of the numerical methods for solving the variational 
equations from the ansatz were increased significantly. The combined effect of these 
three improvements allows us to establish the upperbound on the actionprofile, 
which turns out to be a rather delicate affair. Having established this upperbound 
does not rigorously prove the existence of a new instanton, but makes it very likely 
in our opinion. In that case, we also have, from our results for the non-contractible 
loop, an approximation of the exact solution I*. 

The present paper is primarily concerned with the existence of a new classical 
solution, and we refer the reader to Jl| for the physics that motivates our search. 
The outline of the present article is as follows. In sect. 2 we discuss the constraint 
on the scale and what solutions precisely we are after. In sect. 3 we outline a general 
strategy for the search of new solutions, which is to find non-contractible loops in 
configuration space. One particular non-contractible loop is presented in sect. 4. 
First, we discuss, in subsect. 4.1, the basic ingredients of our construction; then, we 
give, in subsect. 4.2, the ansatz in full detail. The general behaviour of the action 
over this non-contractible loop is discussed in subsect. 4.3 . The actual numerical 
results for the actionprofile are presented in subsect. 5.2. These results show that, 



actionprofile holds, which is the main result of this paper. As a byproduct we obtain 
some numerical results for the well-known BPSTH instanton ||, g]; these results are 
given in subsect. 5.1 . Finally, we return, in sect. 6, to the possible existence of the 
I* solution in SU(2) Yang-Mills-Higgs theory and compare this with the situation 
in other theories. We also discuss, very briefly, potential applications of I* to the 
physics of the electroweak interactions. There are two appendices. In the first of 
these, Appendix A, we give an outline of the calculations behind the numerical 
results of sect. 5. In the second, Appendix B, we describe a similar calculation for 
a static, 3-dimensional solution, the new sphaleron S*, which is directly related to 
the new instanton I*. The reader who is not interested in the details is advised to 
concentrate on sects. 2 and 3, possibly subsect. 4.1, and sect. 6. 



2 Constrained Instantons 



We consider a Yang-Mills-Higgs theory with SU(2) gauge fields W and a single, 
complex Higgs doublet $. The euclidean action of this theory is 



^4ymh = / d 4 x 
Jr 4 



^Tr^ + |D M $| 2 + A(|$| 2 -^y 



where = d„W v - d u W» + [W^ W v ] , = d„ + W„ and W„ = W% a a /{2i) , with 
a a the standard Pauli matrices. Here, and in the following, Greek indices run from 
to 3 and Latin indices from 1 to 3. The masses of the three W vector bosons and 
the single Higgs scalar H are Mw — gv/2 and Mh = \/2Xv. Finite action gauge 
field configurations can be classified by the Pontryagin index (or topological charge), 
which takes on integer values and is defined by 



2 Tr (2 CkA ^ WkX W>iU 



(2) 



A simple scaling argument shows the absence of non-singular 4-dimensional clas- 
sical solutions for the action (|l|) : any configuration can lower its action by collapse to 
a point. For this reason one introduces a constraint on the scale and later integrates 
over the corresponding collective coordinate p. A convenient way to implement the 
constraint is to require || 

J d 4 x O d = 8n 2 c p"- d , (3) 

where = 0,^ (W, $) is a field operator with canonical mass dimension d > 4. 
With a Lagrange multiplier k one then looks for the stationary points of 

^4ymh + k (V" 4 J d 4 x O d - 8 tt 2 c) . (4) 

Concretely, we proceed as follows. First, we have to chose the constraint operator 
Od- The actual choice is irrelevant in the end, as long as we integrate, in the path 
integral, over the collective coordinate p (with the appropriate Jacobian). For purely 
technical reasons, to be explained later on, we have chosen to work with the operator 

8 = [~ Tr 1 e KXilv W KX . (5) 

Having made this choice, we solve the field equations that result from variations 5W 
and 5$ in the constrained action 



for a given value of the dimensionless coupling constant k. Let us denote the fields 
of one particular solution by W* and $*, with a corresponding constrained action 
A*. The scale of this solution is then determined by 



.8tt 2 c 

and its Yang-Mills-Higgs action by 



d 4 x 8 (W* 



1-1/4 



(7) 



A^ MR = A YMH (W\^) . (8) 

These values p* and Ay MR are functions of the coupling constant k, as the fields 
W* and $* are. Eliminating k one obtains ^4ymh as a function of p*. 

This procedure is entirely straightforward. The problem is to discover the solu- 
tions of the field equations from the constrained action (|6|). Any such 4-dimensional 
euclidean solution will be called a "constrained instanton" . [] The prototypical con- 
strained instanton is the solution of Belavin et al. and 't Hooft || |3|. We denote this 
basic instanton by Ibpsth ; or I for short. Its fields are given by (using differential 
forms) 




(9) 

with x 2 = x^X/j,, x M = x^/x and <r M = (l,ia m ). The instanton has Pontryagin index 
Qi = 1. By reflection there is also an anti- instanton I, with equal action but opposite 
topological charge Q\ = —1. Only in the limit p, A — > 0, is it possible to find an 
analytical solution for the radial functions in the ansatz [0, |J 

.2 



/] 



X 



x 2 + p 2 



hi = Jfi, (10) 



with a corresponding action 



1 , . . , 2 \ 87T 2 



^4ymh = ( 1 + 2 (P M w) ) -j + O (A) . (11) 



1 In pure Yang-Mills theory "instantons" are sometimes meant to refer exclusively to solutions 
of the first-order self-duality equations, which then solve the second-order field equations by the 



In sect. 5.1 we will give some numerical results for the action at finite values of p 
and A. 

In this article we look for constrained instantons that are not related to the 
BPST instanton (\Q\ = 1) [Q] or other self-dual solutions (\Q\ > 1) || of pure Yang- 
Mills theory. In fact, our search is for solutions in the vacuum sector (Q = 0) of the 
constrained Yang-Mills-Higgs theory (|B|). No such solutions are known at present. 



3 Strategy 



Explicit construction of instantons in the vacuum sector of the constrained theory is 
not feasible at present. Instead, we present a topological argument for the possible 
existence of at least one such solution. This topological argument involves the con- 
struction of a suitable non-contractible loop (NCL) of 4-dimensional configurations 
of the fields W and <£>, starting and ending at the classical vacuum. In addition, we 
may hope to gain some insight into the structure of this solution, preliminary to an 
explicit construction of it. 

The presence of non-contractible loops in configuration space implies, as Taubes 
M has shown in a somewhat different context, the existence of new solutions of 
the classical field equations. The intuitive idea is that by "shrinking" the NCL it 
gets "stuck" at a point in configuration space, which corresponds (or is close) to a 
stationary point of the action, i.e. a new solution of the field equations. Specifically, 
this is a mini-max procedure, where we take the maximum action on the NCL and 
try to minimize that action. It is essential for this topological argument that the 
euclidean actiondensity @, with positive coupling constants A and k, is a positive 
semi-definite functional of the fields, or, in other words, that it is bounded from 
below. In mathematics this general subject is called Morse theory, which aims to 
relate the critical points of a functional to the topology of the function space on 
which the functional is defined. The method of finding these stationary points by 
a mini-max principle on non-contractible loops goes under the name of Ljusternik- 
Snirelman theory. References to the mathematical literature can be found in ||. 

For the case at hand there is, however, one obvious loophole in the argument. 
It could be, namely, that the mini-max procedure leads to an approximate solution 
consisting of the BPSTH instanton solution I and anti-instanton I at infinite separa- 
tion. In that case there would be no genuine new solution. Clearly, this possibility 
is ruled out if we are able to construct a particular NCL for which 



where oj parametrizes the position along the loop. This upperbound on the action- 
profile is a necessary condition for an existence proof of I* and the main goal of this 
article is to establish it. 




(12) 



4 Non-contractible loop 



In this section we present an ansatz for a non-contractible loop of configurations. We 
start with the basic ingredients that go into the construction, then give the details of 
the ansatz, and, finally, discuss the expected behaviour for the action over the loop. 
The actual profile of the action over the NCL has to be determined numerically, 
these results will be presented in sect. 5. 

4.1 Basic ingredients 

The first step in the construction of the NCL is to give the structure of the fields at 
infinity (|a;| — > oo). This will be based on a non-trivial mapping 

U : S 3 x S 1 ! -> SU{2) ~ S 3 , (13) 

where the first £3 refers to the hypersphere at infinity, Si to the loop of configurations 
and SU(2) to the gauge group. The mapping should belong to the non-trivial 
homotopy class in 774(5*3) = Z 2 , so that, later on, we have a loop of configurations 
that is indeed non-contractible. A specific choice for the SU(2) matrix U is the 
following 

exp[(cu + 71-/2)203] (x ■ cr) exp[— (uj + 71-/2)20-3] (x ■ cry 
Xp/ \x\ 

(l,ia), (14) 

with a loopparameter uj G [— 7r/2, +7t/2]. The fields at infinity (pure gauge, of 
course) are given by 

= -duu- 1 

The second step is to extend these fields inwards. For this purpose we introduce 
two radial functions f(x) and h(x), which approach 1 at infinity and vanish at the 
orgin, in order to ensure continuity. Also, we extend the range of the loop parameter 
uj to [— 7r, +7r] , and make it into a real loop, starting and ending at the same point, i.e. 
the vacuum at uj = ±7r. In this way we arrive at the following NCL of configurations 



U = 

Xfx — 
CTu = 




u G [-7T, -tt/2] U [+7r/2,+7r] : W = 

$ = (h sin 2 + cos 2 10 

uj G [ — 7r/2,+7r/2] : W = -f dU U' 1 

$ = [ C I , (16) 




where the radial functions f(x) and have the boundary conditions 

lim /, h = 1 
/(0) = h(0) = 0. (17) 

The precise form of these radial functions f(x) and h(x) is arbitrary in principle, 
but the mini-max idea is to optimize them at the maximum point (u = 0) on the 
NCL. This is done by solving the variational equations for f(x) and h(x), that result 



from inserting the ansatz (|T^) with uj = in the constrained action (|(J) and making 
variations 5f and Sh. 

The NCL ( pp is quite elegant, but not good enough for our purpose. The reason 
is that the actionprofile A(u) has a maximum value A(0) definitely greater than 2A i. 
Hence, the inequality ( |T2"D does not hold for the simplest possible NCL. This brings 
us to the third, and final, step in the construction of a suitable NCL. We start from 
the observation that there are really two "cores" in (JT^), each of which resembles 
the (anti)instanton (I) I as given in @. The idea now is to separate these cores, in 
order to profit from the attraction of the long-range fields present in the theory, the 
Higgs field in our case. This improved NCL, parametrized by uj G [— 3tt/2, +37r/2], 
has the following structure 

— 37r/2, — 7r] build up of an initial Higgs configuration; 

— 7r, — 7r/2] creation and separation of an II pair; 

— 7r/2, +7t/2] relative isospin rotation of the II pair; 

+7r/2, +7r] collapse and annihilation of the II pair; 

+7T, +37r/2] clean up of the remaining Higgs configuration. 



Ia : 


LU G 


Ila : 


uj G 


III : 


uj G 


lib : 


UJ G 


lb : 


UJ G 



The explicit construction of the NCL is somewhat involved and will given in the 



4.2 Ansatz 

The configurations of the NCL are given by 

I : W = 

$ = (h cos 2 uj + sin 2 uj) 



V2 1 



II, III 



W 



-f dU U- 1 



V2 \ 1 
with the following SU(2) matrix 

U = exp{ttia 3 ] (y_ ■ a) exp[-Qia 3 ) (y + ■ a) 1 " 

and notation 



(19) 



Vtt± 


= y^±/ \y±\ 


o» = 


Cm) 


Vm± 




yo± = 


x ± = x ± D/2 


x\ 


= %l± + r 2 




x 2 m =p 2 + z 2 


z 


= x 3 


t = 


Xq . 



Here D = D{u) determines the core distance and Q = Q(uj) their relative isospin 
rotation 
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(20) 



The functions /, h and g± in ( P| [H5| ) are taken to be axial functions f(r,t), h(r,t) 
and g±(r,t), with the following boundary conditions 



|a;|— >oo 

f(Q,±D/2) = h(0,±D/2) = 



(21) 



and reflection symmetry 



/(r,t) 



f(r, ~t) 

h(r —+\ 



The Pontryagin index (|2]) vanishes for all configurations of the NCL, precisely be- 
cause of this reflection symmetry of the ansatz. This completes the basic construc- 
tion of the NCL. We remark that the only difference compared to our previous paper 



1 1]] is the presence of the functions g± in the matrix U (19), but this change will turn 
out to be essential. 

It remains to specify the four axial functions f(r, t), h(r, t) and g±(r, t) that enter 
the ansatz. Just as for the simple NCL of the previous subsection, the mini-max 
procedure would be to insert the ansatz (|I8| ), for to = and D = d max , into the 
action (|6]) and solve the variational equations for /, h and g±. This turns out to 
be prohibatively difficult and, instead, we make an explicit choice for g±. As will 
be explained in the next subsection, we want the cores to become independent for 
large values of their separation D. This can be achieved by the following choice, for 
example, 

{x%D*/AT + (gg ^ 
9± (x 2 T D 2 /4) a + {(5p i) 4q + {D/2Y a 1 } 

with x\ defined below ( |i~9"D and the parameter values a = 4.0 and (3 = 1.5, obtained 
by trial and error. Here p \ — p i(/c) is the scale for the BPSTH instanton I, see sect. 
5.1 below. Having fixed the functions g±, we can solve numerically the variational 
equations for / and h. We will do this for both uj = and tt/2, and for arbitrary 
values of the distance D. These solutions will be called f(cj, D) and h(u, D), where 
the dependence on the spatial coordinates r and t is implicit. In order to keep the 
variational equations as small as possible, we have chosen the constraint operator 
(||), and not, for example, 

which would give significantly larger expressions. Furthermore, we set the constant 
c in the definition of the scale (0) to the value 

c = 288/21 , (24) 

in order to match the scale of the BPSTH instanton, which will be determined in 
sect. 5.1. 

With the solutions f(cu,D) and h(u,D) in hand, we can, at last, specify the 
functions / and h that enter the ansatz (|i~8|) for the NCL 



I, II : f = f(n/2,D) 
h = h(n/2,D) 



Ill : / = /(0,d max ) cos 2 ^ + /(vr/2,(i max ) sin 2 ^ 

h = h(0, d max ) cos 2 10 + h(n/2, d max ) sin 2 u , (25) 

with the distance D = D(oj) given by (0). Actually, we only need the u = n/2 
solutions / and h down to some small value D = d m in- We can then close the 
loop with a simple deformation of these functions and keep them non-singular, as 
explained in our previous paper JIJ and in Appendix B at the end of this one. By 
choosing d min small enough, the action can be kept arbitrarily small. 

To summarize, the field configurations of the NCL are given by (18 -|20Q with the 



axial functions ( p3| , |2~5|) , where / and h are the solutions of the variational equations 
for uo = or it/2 and D < <i max . 



4.3 Actionprofile 

Our main interest lies in the actionprofile over the NCL. Here, we will give a general 
discussion of the possible behaviour, in order to prepare the way for the numerical 
results to be presented in sect. 5. 

The profile of the constrained action over the NCL is essentially determined by 
the behaviour of the constrained action A(u, D) for the solutions / and h at uj = 
or 7r/2. In fact, / and h solve the variational equations, so that the action attains its 
lowest value, within the ansatz, precisely for these functions. On physical grounds, 
we expect the following behaviour for A. For D — > oo the constrained action should 
approach twice that of the BPSTH instanton I, whereas for D — ► the effects of 
the Yang-Mills interactions should become important, which can be either repulsive 
(u = 0) or attractive (u = vr/2). The Yang-Mills interactions for large values of 
D get suppressed exponentially, with a length scale set by M^ 1 . In this paper, we 
consider the case of vanishing Higgs mass = or A = 0. This means that at 
very large distances D, where the Yang-Mills interactions drop out, only the effect 
from the Higgs field remains. Moreover, the Higgs interaction at large distances is 
attractive and basically independent of the relative isospin rotation Q of the cores. 
So, the simplest behaviour we expect is the one sketched in fig. la. At this moment, 
we can also explain the need for the additional functions g± in the ansatz (|19|) . For 
g± — 1, namely, there would be "tidal effects" from one core on the other, which 
decrease rather slowly with the distance D. If, instead, g + vanishes approximately 
near the t = +D/2, and g_ near t = —D/2, the cores become independent faster 
with increasing distance. This. then, is the reason behind our choice fEHT) for these 



As mentioned before, the actionprofile over the NCL is essentially determined by 
the behaviour of A(0,D) and A(n/2, D). This follows from, in particular, (p0|) and 
( |25|) above. Concretely, the actionprofile, which is an even function of u, is obtained 
as follows. We start in fig. la on the u = tt/2 curve at D = and move out to 
D = d max , then we go straight up to the to — curve and back again, and finally 
we return in the same way to D — 0. We have verified that the action over part 
III of the NCL, with the functions (^5|), has indeed a single maximum at ou — 0. 
From the curves of fig. la, then, we would conclude that the maximum action over 
the NCL can remain below 2A i, provided we have a large enough value of c? max . In 
fact, the optimal choice (i.e. lowest possible maximum action) would be to move 
out to a distance d* with a corresponding maximum action A*, both of which are 
indicated in fig. la. However, another possible behaviour of A(0, D) and A(^/2, D) 
is sketched in fig. lb. In that case, the maximum action would always stay above 
2A i, regardless of the value of (i max - A priori, there is no way to decide between the 
two types of behaviour shown in fig. 1 and we need an explicit calculation to settle 
the matter. 



5 Numerical results 



We present here some numerical results for the NCL constructed in the previous 
section. These results are for the case of vanishing quartic Higgs coupling constant 
A = 0. We also give some results for the BPSTH instanton I. Distances will be 
expressed in units of M^ 1 and the action in units of 8n 2 /g 2 , which is the action 
of the BPST instanton in pure Yang-Mills theory. A brief outline of the numerical 
methods is given in Appendix A. 

5.1 BPSTH instanton 

We have solved numerically the variational equations for the radial functions / i(x) 
and h i{x), that result from inserting the BPSTH ansatz (^) into the constrained 
action (^) and making variations 5f i and 8h j. As explained in sect. 2, these 
solutions / i and h i depend on the coupling constant « of the constrained action A. 
We determine the scale p of the instanton from (^), with the constant c = 144/21, 
and the Yang-Mills-Higgs action Aymr from (R), where W* and $* now refer to the 
ansatz @ with the functions / 1 and h i. This particular choice for the numerical 
value of the constant c in (^) reproduces the scale p that enters the analytical 
solution fllOD of pure Yang-Mills theory. Our numerical results are collected in Table 
1, where, for future reference, we also give results for X/g 2 = 1/8 (M H = M w ) and 
X/g 2 = 100/8 (Mh = 10 M w ). The A = results for A Y mh vs. p are compared, in 
fig. 2, with the expression fli~T]) , which is valid in the limit p — > 0. For values p ^ M^ 1 
the radial functions of the numerical solution differ significantly from the analytical 
functions (|10D, resulting in a lower value of the action. Another consequence is that, 
for finite values of the scale p, the gauge fields are no longer self-dual. 

The numerical results for the instanton I are relatively easy to obtain, as the 
variational equations for f i and hi are ordinary differential equations (ODEs). 
These results will serve as a check on those of the NCL, to which we turn now. 

5.2 Non-contractible loop 

We have solved numerically the variational equations for the axial functions f(r,t) 
and h(r, t), that result from inserting the ansatz of the NCL (|i~8f) , for uj = or tt/2, 
into the constrained action (||) and making variations 8f and Sh. These variational 
equations consist of two coupled, non-linear partial differential equations (PDEs), 



solutions / and h, we also determine the Yang-Mills-Higgs action Aymh and the scale 
p from the expressions ([!]) and (0), respectively, with the numerical value ( p4|) for 
the constant c and the NCL configurations for W* and $*. As discussed in subsect. 
4.3, we are primarily interested in the behaviour of the constrained action of the 
solution as a function of the core distance D. For D — > oo we reproduce the results 
of the previous subsection for 2 A j (Table 1). In fig. 3 we give the numerical results 
for the constrained action, normalized to its asymptotic value, for a relatively large 
value of the constraint coupling constant, namely k — 1. Figure 4 shows the same 
results on an expanded scale. These numerical results seem to agree, for k = 1, 
with the simple behaviour sketched in fig. la. In contrast, the results for k = 10~ 3 , 
shown in fig. 5, display the alternative behaviour of fig. lb. This different behaviour 
for small and large values of the scale was also seen in our previous results []IJ , but 
we cannot directly compare the scale p there with the scale p here. 

We are reasonably confident that the results of figs. 3-5 are reliable, for the 
following reasons. First, the results are stable when the calculation is repeated 
on different grids (ranging from 25 x 50 to 100 x 200 points) and with different 
compactified coordinates for r and t, see Appendix A. The error of the data points 
in fig. 4 is estimated to be less than 0.1%, which is about a factor of 10 better than 
the numerical accuracy of our previous results U . Second, the variational equations 
for / and h are solved by relaxation, which implies that the exact solution can have 
a somewhat lower action. In other words, the "dip" of the to = curve in fig. 4 can 
only get deeper. Third, and as is often the case with large numerical calculations, 
we beleive our results because they behave in the way we expect them to do. In 
particular, the dip of fig. 4 becomes more shallow with X/g 2 increasing. Also, we 
have an heuristic understanding of the origin, for small values of k, of the "bump" 
on the oo = n/2 curve in fig. 5. Fourth, a similar calculation for static fields gives 
comparable results, see Appendix B. 

To summarize, we have constructed a NCL, for which the important upperbound 
( |T2"D holds, provided the scale p is large enough. Specifically, we obtained for X/g 2 = 
and k — 1 (see figs. 3-4) 

maxv4 NCL H = 1-994 A 1 < 2 A i . (26) 

The optimal maximum (u = 0) configuration of the NCL gives, moreover, an ap- 
proximation of the conjectured constrained instanton I*. This configuration leads 
to the following estimates for I* 



cf ~ 10 

^ymh ~ 4.2 8n 2 /g 2 . (27) 

In fig. 6 we show the corresponding actiondensity Oymh^, t) and Pontryagin-density 
q(r,t), averaged over the polar angle 9, see Appendix A. This shows that our con- 
figuration is still a very loose molecule and we expect the exact solution I* to be 
tighter and more cigar-like perhaps. But it is also clear, from fig. 4 especially, that 
the Yang-Mills cores are very hard and that d*, which is in essence the distance 
between the points of vanishing Higgs field, cannot be much smaller than the width 
(~ 2 p*) of the configuration. 



6 Discussion 



We have constructed in this paper a non-contractible loop (NCL) with a maximum 
constrained action less than twice that of the BPSTH instanton, provided the length 
scale is fixed at a large enough value. This was established for the case of vanishing 
Higgs mass M H = 0. We expect it to be possible to extend the result to all values 
Mn < My/. Note that in the full electroweak theory there is also the photon field, 
which can provide the necessary attraction if the Higgs field becomes too short of 
range. As it stands, this upperbound on the actionprofile over the NCL is only a 
necessary condition in an eventual existence proof. Still, we are optimistic about the 
existence of I* in SU(2) Yang-Mills-Higgs theory. There exist, of course, analogous 
static solutions (sphalerons) in Yang-Mills-Higgs theory || [7], ||, but there are also 
encouraging results on instantons in pure Yang-Mills theory, which we will now 
discuss. 

There has been a long-standing conjecture [||, based on the analogy with har- 
monic maps from 5*2 to S2, that all solutions in euclidean SU(2) Yang-Mills the- 
ory over 5*4 are necessarily self-dual or anti-self-dual (W^ u = i^e^xWux). In the 
vacuum sector, in particular, there would be no other solutions besides the classical 
vacuum itself. This conjecture has been proven false recently. In fact, Sibner et al. 
[|IIJ showed the existence of infinitely many non-self-dual solutions in the vacuum 
sector, with an action 

A YM = m 16 7T 2 /g 2 + AA YM , (28) 
for integers m > 2 and AAym < 0, which, most likely, depends on m also. Later, 



several solutions were constructed explicitely by Sadun and Segert The exis- 



tence proof of Sibner et al. [T(| goes by a mini-max procedure over non-contractible 
loops, where an important ingredient is the inequality contained in (|28|), which is 
analogous to ours flT2"|). The other main ingredient is an equivariant weak com- 
pactness theorem, where equivariant refers to an U(l) symmetry of their ansatz. 
The integer m, which appears in (|28|), labels the embedding in SU(2) of this U(l) 
symmetry (rotation angle 7 G [0, 27r]), namely by the matrix exp[m7<J3/(2i)] for the 
gauge field transformations. 

Equivariance is known to be a powerful tool in Morse theory. For this reason, it 
may be of importance to note that our ansatz ([18|, pT9| ) also has an U(l) symmetry, 
and precisely for the case m = 1 excluded in the existence proof |H| and the explicit 
construction [TT| of the pure Yang-Mills solutions. The difference is that our theory 



for the crucial upperbound ( O ) on the action. The constrained action (^j) of this 
theory now has three terms extra, compared to the case of pure Yang-Mills theory. 
The constraint term, in particular, is designed to prevent collapse and, naively, we 
see no way how a regular solution I*, related to the presence of non-contractible 
loops in configurations space, could fail to exist. 

Just as the sphaleron S [0, is associated with the BPSTH instanton I (loosely 
speaking S is a constant time slice of I), we expect a new sphaleron S* [0 to 
be associated with I*. In Appendix B we give some numerical results for a non- 
contractible sphere of 3-dimensional configurations, which support the existence of 
this new sphaleron S*. After these results were obtained, we succeeded in construct- 
ing the solution S*, on which we will report elsewhere. We expect the construction 
of I* to proceed in the same way, only with greater technical complications. Hence- 
forth, we take for granted the existence of the new constrained instanton I* in SU(2) 
Yang-Mills-Higgs theory. 

An interesting question is to see what happens when massless fermions are intro- 
duced into the theory, or possibly fermions with Yukawa couplings to the Higgs. We 
then expect I* (and S*) to have fermion zero-modes. The reason is that, as explained 
in our previous paper jl], there is spectral flow of the Dirac eigenvalues along the 
NCL fllSp. This spectral flow was indirectly monitored in the numerical calculations 
of sect. 5 and our results support the claim that I* has fermion zero-modes. Conse- 
quently, there should be a new (B + L conserving) effective fermion vertex from I*, 
with double the number of lines compared to the one from the BPSTH instanton 
[[|. We will now turn to different instanton solutions, which may be related to I* in 
one way or another. 

The pure SU{2) Yang-Mills solutions fL0|, [D] discussed above can be expected to 
have counterparts in SU(2) Yang-Mills-Higgs theory, i.e. there will be corresponding 
constrained instantons. In general, these solutions will have a rather large action, 
since for them we have already that ^4ym ~ 32 ir 2 /g 2 , whereas Aym ~ 16 n 2 /g 2 for 
I*. We expect I* to be the lowest action non-trivial solution in the vacuum sector 
of the theory. More interesting, perhaps, is the inverse question, wether or not 
the constrained instanton I* of SU(2) Yang-Mills-Higgs theory has a corresponding 
solution in pure Yang-Mills theory. We conjecture there to be no such counterpart 
in pure Yang-Mills theory, because we see no obvious source of attraction. The fact 
that we found the upperbound on the action (|l^) to be violated for small values of 
the constraint coupling constant (see fig. 5) is suggestive, but is not really conclusive, 



Different Yang-Mills-Higgs theories can also be considered. If W* and $* are 
the fields of the I* solution in SU (2) Yang-Mills-Higgs theory, then it is possible to 
embed them into a larger theory, provided that SU(2) is a subgroup of the gauge 
group G and that the SU{2) Higgs doublet (together with its vacuum expectation 
value) can be embedded in the larger Higgs representation. An example is £77(3) 
Yang-Mills-Higgs theory with a complex triplet of Higgs. The embedding is then 
given by 
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(29) 



where the crosses indicate, symbolically, the SU (2) solution. Note that these embed- 
dings can lead to unexpected solutions of the field equations, since the larger theory 
may not even have non-contractible loops (it 4(G) = 1, as is the case for G = S77(3), 
for example). The NCL was a tool to find the S77(2) solution, but once we have 
found the solution, we can forget about the NCL and simply verify the fact that the 
ansatz solves the field equations, which is then carried over to the larger theory. 

Finally, we comment on possible applications of the conjectured new instanton 
solution I*. Of course, S77(2) Yang-Mills-Higgs theory is at the heart of the Glashow- 
Weinberg-Salam model for the electroweak interactions. Evaluating euclidean path 
integrals for electroweak processes, it may be important to know that in addition to 
the classical vacuum there is a new stationary point, i.e. the constrained instanton 



I*. It is well-known |U| [L4]] that such stationary points (and vacuum instability, in 
general) can play a role in the asymptotics of perturbation theory. Furthermore, 
this new stationary point could contribute directly to the euclidean path integrals of 
certain forward elastic scattering amplitudes, which control the total cross-sections 
for the corresponding processes. In our previous paper |l| we have shown that 
this contribution, evaluated classically f], suddenly becomes important as the parton 
center-of-mass energy y/s increases. In fact, this threshold energy is determined by 
the structure of the I* solution, and from our approximation of that solution we 
obtain 



threshold 
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^YMH / ^YMH 
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d* M w 



(30) 



2 For a semiclassical calculation we also need to know the negative modes around the classical 
solution. In pure SU (2) Yang-Mills theory it has been shown Jl5| that there must be at least two 
negative modes, instead of one; wether or not the same holds in Yang-Mills-Higgs theory is not 
known at the moment. Note that even a new, "non-perturbative" contribution to the real part 



where we have used the definition 



which is close to the true value 3.04 M w /a w for the sphaleron energy at A = 0. With 
the numerical estimates given in (p7|) we find, not unexpectedly perhaps, that the 
threshold in the parton center-of-mass energy is of the order of E$ ~ 10 TeV. We 
intend to discuss the applications of I* in a separate publication. 
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Appendix A : Numerical methods 



We decribe in this appendix the algebraic and numerical calculations for the results 
reported in sect. 5. These calculations are straightforward, but cumbersome. Hence, 
we will give the main points only and leave many technical details out. In the first 
part of this appendix, we review the algebraic calculation of the actiondensity for 
our ansatz. In the second part, we discuss the numerical solution of the variational 
equations from this actiondensity. 

Al : Algebraic calculation 

The main algebraic calculation consists of two steps. The first step is to insert the 
ansatz for the NCL into the constrained action. We use axial coordinates 

x o = t , 

£3 = z = r cos 9 , 

X2 = p cos (j) = r sin 9 cos <fi , 

Xi = p sin0 = r sin 9 sin0 , (A.l) 

and make all distances dimensionless with M^ 1 . The total constrained action (|6]) 
for the ansatz (|T8|) takes the form 

I /*oo />oo r~ r2n 

A=— dt drr 2 d9sin9 d(f> a , (A.2) 
g 2 J -00 Jo Jo Jo 

with a rotationally invariant actiondensity 

a = a(z, p, t) . (A. 3) 



As a matter of fact, it is already clear from the ansatz (18 , |19| ) that a rotation of 
(xi, X2) can be compensated by a global gauge transformation. 

The second step is to perform two of the integrals in the action ( |A.2j ). The one 
over the azimuthal angle <fi is trivial, of course. The integral over the polar angle 9 
can also be performed, since the functions /, h and g± depend, by construction, on 
r and t only. In addition, we have the reflection symmetry t — ► — t, so that the final 
expression of the constrained action for our ansatz takes the form 



fOO /*00 



A = — — / dt dr — d. (A.4) 



<T JO J 7T 

The averaged actiondensity a = a(r,t) has the following structure 



(drfY + Zymu (dtfd r f) + Zymoo (/(l — /)) , 



a H = 2(^) 2 + 2(9 r /i) 2 + ^Hoo(Ml-/)) 2 + 4-(/ i 2 -l) , 

a c = K(Zc2o(d t f) 2 + Z c0 2(d r f) 2 + Z cu (d t fd r f))(f(l-f))\ (A.5) 

where d r and <9i denote partial derivatives with respect to r and t. The eight co- 
efficients Z are complicated rational functions of the variables r and t, together 
with the functions g±{r,t) and their various partial derivatives. These coefficients 
depend also on the loop parameter u. With some further effort, ( |A.5| ) can be put in 
manifestly positive definite form. We have used the symbolic manipulation program 
FORM Jl6| for the algebraic calculation of the actiondensity (|A.5| ). 



A2 : Numerical calculation 



We now have to solve numerically the variational equations from the action ( |A.4| ). 



We proceed in three steps. The first step is to compactify the coordinates r and 
t to the variables x and y, respectively. Since we are interested in the long-range 
behaviour of the Higgs field, we choose for y{t) a rather slow dependence on t, 
specifically 

t-tr 

y = y c + (i - y c 



l + t - 1 



c i 



y c = t c /(l + 2t c ) 

t EE t/t scale (A.6) 

and similarly for x(r). Furthermore, we let t c correspond to the core position t c = 
D/2 and set, typically, r c = pi and r scale = t scale = 2 p h where p x = p i (k) is the 
scale of the instanton I, see sect. 5.1. It is straightforward to make these changes of 



variables in the action (|A.4j ), and we write the result as 



8vr 2 r 1 r 1 

A = — — / dy dx a , (A. 7) 

g 2 Jo Jo 

with the Jacobians absorbed into a = d(x, y). We look for the two functions f(x, y) 
and h(x, y) that minimize this action, with mixed Dirichlet and Neumann boundary 
conditions as indicated in fig. 7. 



The second step is to discretize the integral ( |A.7| ). We use a rectangular grid for 
x and y 

x = iAx i = 0, . . . , I 

y = jAy j = 0,...,J (A.8) 

and write for the functions at the lattice points 

f(iAx, jAy) — > fij . (A.9) 

Furthermore, we use central differences for the partial derivatives of / and h, for 
example 

OJVAr.jAy) -> : • (A.10) 

The resulting discretized action is, however, numerically unstable : minimization 
leads to functions fc j and h it j that take on alternating values of approximately 
and 1. We have chosen to employ the following two countermeasures. First, we 
"double" the boundary conditions at infinity, namely 

hj = hi-ij = 1 , (A. 11) 

and similarly for /. Second, we "smear" the constraint term dc{x,y) in the action- 
density Specifically, we take for the first term in a c , cf. (|A.5|), 



kZc20 { 2Ay ) { 2 ) ^ 2 ) ' (A ' 12) 

and similarly for the other terms. In this way we end up with a discretized action- 
density dij at the gridpoint For the total constrained action we have 

8vr 2 1 J 

A = — zl Ax Ay Wij a iyj , (A. 13) 

9 i=o j=o 

with Wij the weightfactors of, for example, the extended trapezoidal rule. 

The third, and final, step is to solve numerically the variational equations for 
fij and hi j coming from the discretized action (|A.13| ). These equations are highly 
non-linear, especially for the case k = 1 we are interested in, and we need a method 
that can handle this. We have succesfully employed the method of non-linear over- 
relaxation (NLOR) 0, using grids of, typically, 25 x 50 points. Our FORTRAN 
program starts with some trial functions for /y and hij, which are then relaxed, 
first with NLOR, but in the end also with some sweeps of underrelaxation, if neces- 
sary. In this way we obtain smooth configurations / and h, with a definite value for 

+- Vi o nnn of roinorl o r>4~ i nn T 1 !-! o ovo r>4" \1 nf inn m Qir a~t ran It o ~i ro o cnm our In o ^ 1 mirov o i nn 



Appendix B : Non-contractible sphere 



We consider in this appendix static, 3- dimensional configurations of the gauge field 
W and the Higgs field <3>. The energy functional for these fields is 



-Eymh = / d 3 x 



•^Tr^ + |An$| 2 + A(W-y) 2 



;b.i) 



with the same notation as in (|l|). Henceforth, we set the quartic Higgs coupling 
constant A = and, for brevity, refer to -Eymh as E. 

It is well-known, by now, that a non-contractible loop (NCL) of static configura- 
tions leads to the existence of a static, but unstable, classical solution, the sphaleron 
S [[/], |§. It is not difficult to construct also a non-contractible sphere (NCS) of static 
configurations ||12|| . A mini-max procedure over this NCS suggests the possible ex- 
istence of a new sphaleron S*, provided we exclude the case of two sphalerons S at 
infinite separation. This loophole is closed if we are able to construct a NCS, for 
which 

max £ N cs(/i, v) < 2 E s , (B.2) 

where \i and v parametrize the position on the sphere. Evidently, this discussion 
parallels the one for the new instanton I* in the main part of this paper. Moreover, 
the ansatz of sect. 4 can easily be adapted to the present case. We will simply state 
the resulting ansatz for the NCS and give our numerical results, which establish the 
important inequality (|B.2|). These explicit results support the somewhat heuristic 



arguments given in our previous paper 12 



The NCS is parametrized by the square v G [— ir, with the boundary ■ 
7T or \u\ = 7r corresponding to the classical vacuum. Writing [fiu] = max(|/i|, \u\ 
the configurations of the NCS are 



tt/2 < \jmv) < 7T : W = 







* = (i-(i-/0HH) 7i \ , 



< [H < tt/2 : W=-f dU U~ x 



^ = h J_ul I (B.3) 



y/2 V 1 

with the following 577(2) matrix U for /z, u G [— 7r/2, +7t/2] 



U = exp[(z/ + n/2)ia^\ exp[(/x + n/2)y^ ■ ia] exp[— (v + tt/2)z(T3] ■ 

exp[-(/i + 7r/2)y + -za] (B.4) 



and notation 



x m± = (xi, x 2 , x 3 ± D/2) = (psin0, pcoscp, z ± D/2) 

ym± = (S'iPsin^, g±pcos(f), z ± D/2) . (B.5) 



The core distance D = D[y) is given by 



vr/2 < |z/| <tt : D = 

0< < 7r/2 : D = <i max cos 5 ^ 



(B.6) 



with 5 > a free parameter. The axial functions f(p,z), h(p,z) and g±{p,z) have 
the following boundary conditions and reflection symmetry 

>oo 

/(0,±£>/2) = /i(0,±D/2)=0 

h(p,z) = h(p,-z) 
9±(p,z) = 9t(Pi~ z ) ■ ( B - 7 ) 

F° r 9±{Pi z ) we take again the functions (£2^), with the same coefficients a and /?, 
and where x± is now defined by (|B.5|) . It remains to specify the two axial functions / 
and h. We will give two alternative constructions, both of which lead to the desired 
inequality (|B.2| ). 



The first construction for the functions / and h in the ansatz (|B.3|) parallels 



the procedure followed for the new instanton in sect. 4, with the loop parameter uj 
there corresponding to the sphere parameter v here. The procedure is to solve the 
variational equations for / and h at fi — and v = or tt/2, for arbitrary values 
of the core distance D. We denote these solutions by /(0, 0, D) and f(0, 7r/2,D), 
and similarly for /i, where the dependence on the spatial coordinates p and z is 
implicit. In fig. 8 we show the corresponding energy values. For v = there is a 
clear minimum at a core distance d* ~ 7 M^} 1 . With these solutions / and /i we 
3 The dip for the instanton case (fig. 4) is more shallow, one reason being the fact that the 



can specify the functions / and h in the ansatz ( |B.3| ) 



0<|H<tt/2 : / = cos 2 v /(0, 0, D) + sin 2 u /(0, 7r/2, D) 
/i = cos 2 1/ fc(0, 0, D) + sin 2 v h(0, tt/2, D) 

-D = rfmax COS 5 1/ 

tt/2 < \u\ < 7T : h = h(0,n/2, 0) . (B.8) 

Choosing <i max = d* and the parameter 5 sufficiently small, we have verified that the 
energy stays everywhere below 2E$. In particular, the energy over slices of the NCS 
at constant values of v have their maximum at p, — and drop to zero monotonically 
for p, — > ±7r. The energy surface, however, is very steep for z/ ~ ±7r/2, and we prefer 
to show it for a somewhat different construction. 

This second, alternative construction for the functions / and h in the ansatz 
( |B.3| ) is as follows. The procedure is to solve the variational equations for / and h 
at p = for the whole range of values v G [— 7r/2, +7r/2], varying the core distance 
simultaneously .D = rf max cos 5 v. Actually, we do not need the complete range 
for v, only the interval [— ^Wn, + z> WaL so that D runs from <i max to some smaller 
value d min (= d max cos 5 z/ min ). We denote these solutions by f(0,v,D;p,z) and 
h(0,u, D; p, z), where now the dependence on the spatial coordinates p and z is 
explicit. For v G [—n/2, — i/ min ] U [+f mm , + 7r /2] we use a simple deformation of the 
functions / and h at v = f mm , in order to move the zeros of / and h appropriately. 
We denote the three segements of the NCS by 



III 

II 

I 



o < M < u, 



mm 



*Vin < \v\ < TT/2 
7r/2< < 7T 



and take for the functions / and h in the ansatz ( |B.3| ) 

III : / = 7(0, v, D; p, z) 
h = h(0, v>, D; p, z) 
D = d max cos 5 v 

p,z + d min /2-D/2) + x 2 _f(0 

j ^min; ^miru 

12 + D/2) 

x+ + x 2 

^min) ^mirij P,z + d min /2 - D/2) + x_ h{0 j ^min) Gemini p, 2 CJ-min 

/2 + D/2) 

X-\- ~\~ X— 



D = d m ax COS 6 V 



d\ava ^(^min) 



j • fi ^(^' ^ mm ' ^minj Pi % ^min/2) "I" h(0, l^raim ^mim Pi % <^min/2) 

This construction leads to a smooth energy surface over the NCS, shown in fig. 
9 for one particular set of parameters. The energy surface over the NCS is, by 
construction, invariant under (p, v) — > (— //, — v). Again, the maximum energy over 
the NCS is reached at p — v — and its value ( 1.94 E$, for the parameters chosen 
in fig. 9) obeys the inequality ( |B.3| ). 

To summarize, we have constructed in SU(2) Yang-Mills-Higgs theory non- 
contractible spheres of static configurations, with energies everywhere below 2 E$. 
This suggests the existence of a new sphaleron S*. Furthermore, we have, for the 
case of vanishing Higgs mass, an approximation of that solution from the optimal 
maximum configuration on the NCS, with a core distance and energy given by (see 
fig. 8) 

E* ~ 1.92 E s , (B.10) 

where E$ takes the numerical value 3.04 Mw/a w . In fig. 10 we show the energy 
density of this configuration. Finally, we note that the 3-dimensional configuration 
( |B.3| ) at /i = v = is essentially equivalent to the slice x% = of the instanton con- 



figuration (|I"8| , [I~9| ) at oj = 0. It is in this sense that the new sphaleron S* corresponds 



to a constant time slice of the new instanton I*. 
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Table 1: Numerical results for the BPSTH instanton (0) in the constrained Yang- 
Mills-Higgs theory (|6|). The classical theory has two dimensionless coupling con- 
stants k and X/g 2 , which control, respectively, the strength of the constraint term 
in the action and the mass ratio of the Higgs scalar H and vector bosons W . The 
scale p of the instanton is given in units of M^ 1 , the constrained action A and the 
standard Yang-Mills-Higgs action Aymr in units of 8ir 2 /g 2 . 



Figure captions 



Fig. 1 : (a) Sketch of the constrained action A as a function of the core distance 
D, for the case of repulsive (u = 0) or attractive (u = n/2) Yang-Mills interactions. 
The Higgs mass vanishes and for large distances there is attraction between the 
cores. The horizontal dashed line corresponds to twice the action of the BPSTH 
instanton. Also indicated are the minimum value A* on the uo = curve and the 
corresponding distance d*. 
(b) Alternative behaviour of A vs. D. 

Fig. 2 : Yang-Mills-Higgs action Aymh (in units of 8TT 2 /g 2 ) of the numerical so- 
lution of the variational equations for the BPSTH ansatz @, as a function of the 
scale p (in units of M^} 1 ). The quartic Higgs coupling constant A vanishes. The 
dashed curve represents the analytical result (|11|), valid for small values of p. 

Fig. 3 : Constrained action A (normalized to its asymptotic value) of the numerical 
solution of the variational equations for the configurations of the non-contractible 
loop (0), as a function of the core distance D (in units of Myj} 1 ). Closed and open 
symbols correspond to, respectively, repulsive (u = 0) and attractive (u = tt/2) 
Yang-Mills interactions. The coupling constants A and k in the constrained action 
(H) have the values and 1, respectively. 

Fig. 4 : Same as fig. 3, but with an expanded scale for the constrained action A. 
Fig. 5 : Same as fig. 3, but for a smaller value of the constraint coupling constant 

K. 

Fig. 6 : Actiondensity oymh(^, t) and Pontryagin-density q(r,t) for the optimal 
maximum configuration of the non-contractible loop (u = 0, D = 10 , A = 
0, K — 1), see fig. 4. Both densities have an arbitrary normalization of their 
maximum to 100. The coordinates r and t are in units of M^} 1 . The complete 
configuration is obtained by reflection t —>■ —t , under which the actiondensity is 
invariant, but the Pontryagin-density changes sign. 

Fig. 7 : Boundary conditions for the functions f(x,y) and h(x,y), defined over the 



and 1 at spatial infinity (x = 1 or y = 1). Neumann boundary conditions, for x = 
or y = 0, are indicated by N. 



Fig. 8 : Energy E (normalized to its asymptotic value) of the numerical solution of 
the variational equations for the configurations of the non-contractible sphere ( |B.3| ), 
as a function of the core distance D (in units of M^ 1 ). Closed and open symbols 
correspond to different parameters of the ansatz, respectively (// = 0, v — 0) and 
(/i = 0, v — 7r/2). The quartic Higgs coupling constant A vanishes. 



Fig. 9 : Energy E (with arbitrary normalization) over the non-contractible sphere 
(B.3, |B.9| ). The energy surface over the whole sphere is obtained by reflection sym- 
metry (/i, v) — > (— //, —i/). The parameters for these results are <i max = 10 Myj} 1 , 
t'min = 77r/16 and 5 = 1. 



Fig. 10: Energydensity e(p, 2) (with arbitrary normalization) for the optimal max- 
imum configuration of the non-contractible sphere (fi = v = 0, D = 7 M^ 1 , A = 0), 
see fig. 8. The coordinates p and z are in units of M^ 1 . The complete configuration 
is obtained by reflection z — > —z. 



